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We develop a non-adiabatic perturbation theory for oscillations involving an arbitrary 
number of neutrino and antineutrino species, including the possibility of lepton number vio- 

■ lation which we treat as a small effect. We interpret the physics of such an approach for the 
C*") | one generation case by introducing the notion of adiabaticity for neutrino and antineutrino 

■ oscillations in analogy to flavor oscillations. We find that in a CP-odd matter environment a 
small lepton number violation in vacuo can be enhanced. Eventually, we apply the perturba- 
tion theory to the two generation case and work out an example for manifestations of lepton 
number violation, which can be solved both perturbatively as well as analytically thereby 
further clarifying the nature of the perturbation expansion. 

Over the last decade or so the physics of neutrino flavor oscillations has seen considerable 
developments. Experiments have found that the three neutrinos discovered so far undergo flavor 
^ ■ oscillations. Two of the three associated flavor mixing angles have been measured, the third one 

has been constrained to be small, but potentially non-zero. Moreover, the two mass splittings 
between the neutrino mass eigenstates have been determined. Despite all progress, a few neutrino 
properties remain elusive. It is, for instance, yet unknown whether there is any CP violation in 
the lepton sector or whether neutrino mass eigenstates obey a normal or inverted hierarchy pQ. 

There are, however, also some experimental findings in neutrino oscillation experiments 
which cannot be reconciled with the standard description of neutrino flavor oscillations. The 
latter findings are therefore often referred to as neutrino oscillation anomalies. Such findings as 
reported by the LSND [2] and MiniBooNE O H] collaborations have given rise to a plethora of 
theoretical speculations how to explain them. The latter range from additional sterile neutrino 
species, non-standard neutrino interactions [5] and extra spatial dimensions [6] to CPT violation 
in the lepton sector [7J . 

Neutrino oscillations depend on neutrino masses. One of the most important questions 
connected with neutrino masses is the possibility of lepton number violation and consequent 
Majorana nature of the neutrinos. Lepton number violation should therefore affect neutrino 
oscillations, and in fact the effects can come in two forms. First of all, in addition to flavor 
oscillations, one can have lepton number violating oscillations, i.e., processes which can be 
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termed neutrino-antineutrino oscillations [8]. Second, the usual flavor oscillation probabilities 
might be modified by the presence of lepton number violating terms in the Lagrangian. 

Despite the fact that the reported neutrino oscillation anomalies might not last and may 
eventually be refuted by future neutrino oscillation experiments, the question of whether physics 
beyond the established picture of neutrino flavor oscillations with similar signatures exists is ul- 
timately an experimental one. In this letter we propose a generic parameterization for lepton 
number violating neutrino flavor oscillations in order to qualitatively understand their ramifica- 
tions and identify deviations from the signals predicted by lepton number conserving oscillations. 

Let us, for such purposes, start by writing down the most general form of the Schrodinger 
equation involving active neutrino and antineutrino fields in flavor space. We denote the left- 
chiral, active neutrino fields collectively by the boldfaced symbol v = (f e , v^, v T , . . .). The active 
antineutrino fields are right-chiral, and they are denoted by v = (p e , u T , . . .). The Schrodinger 
equation has the form 



dt\P(t)J \v(t) 

where H is the Hamiltonian. We assume that the energy of the neutrino beams are such that 
the sterile states, which are expected to be heavy if a see-saw mechanism is assumed to be 
responsible for neutrino mass generation, are not excited and thus dynamically decouple from 
neutrino flavor oscillations. 

In the block form that has been used to write this equation, we can write the total Hamil- 
tonian in the form 

where we have separated out the block-diagonal and off-diagonal parts. Since the fields in v 
are assigned lepton number equal to +1, whereas those in v have the opposite lepton number, 
it implies that Mo is the lepton number conserving part of the Hamiltonian, whereas <5H is 
lepton number violating. While writing the explicit form for SM, we have used the assumption 
that all neutrinos are stable particles, so that the Hamiltonian is Hermitian. This assumption 
also implies that the block matrices H and H are both Hermitian. Note, moreover, that H is 
generically different from H, since neutrino and antineutrino fields can be subject to CP non- 
conserving interactions as is the case, for instance, in elastic forward scattering of neutrinos off 
background fields [9]. 

If CPT is a symmetry of the Hamiltonian, there are additional constraints on the elements 
of M. These constraints can be identified by recognizing that the active antineutrino states are 
CPT transforms of the active neutrino states. In general, if an operator &{t) commutes with 
CPT denoted by the operator 0, we can write 

(0o| 0(t) [06) = {b\ ] {t) \a) = (a\ 0(t) \b)* (3) 

for arbitrary state vectors \a) and \b), and denoting the state Q\b) by \®b). For the active 
neutral fermions which appear in Eq. ([I]), we can choose the ordering as well as the phases of 
the states in such a way that 

e k) = \u a ) . (4) 

Therefore, Eq. ([3]), along with the Hermiticity of the Hamiltonian, would then imply the relation 
(y a \ H(t) \v b ) = {Qu b \ H(t) \Qu a ) = {u b \ W(t) \u a ) . (5) 
In terms of the block-diagonal notation introduced in Eq. ([2]), this can be written as 

H(t) = H T (t). (6) 
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Going through similar steps, we can also show that 

B(t) = B T (t), (7) 

i.e., the block B is symmetric. Of course lepton number violating neutrino-antineutrino oscilla- 
tions can occur even in scenarios with CPT violation, in which case Eqs. © and (J7J) need not 
hold; but also for the case in which Eqs. ([6]) and ([7]) hold, i.e., CPT is an exact symmetry, still 
lepton number violating neutrino oscillations are possible. 

For the purposes that concern us in this paper, it is more convenient to introduce the time 
evolution operator U(t, to) via 

The time evolution equation thus takes the form 

^u(Mo) = -m(t)v(t,t ). (9) 

Henceforth, we will often take to = and will omit it in the notation. Let us have a look at the 
CPT properties of the time evolution operator. CPT symmetry would imply that the evolution 
operator satisfies the condition 

eu(t) = U(-i)0, (10) 
and therefore the analog of Eq. ([3]) would be 

(@a\ U(t) |66) = (o| U(-i) \b)* . (11) 

In particular, it says that 

faM)) = WbWa(-t))* = WaWb{t)) , (12) 

implying that the oscillation probability of an v a going to is the same as that of a going 
to v a : 

P{v a ^v h -t)=P{u b ^v a -t). (13) 

In addition, Eq. (fTTj) also implies that 

P(y a ^v h -t) = P{y h ^v a \t). (14) 

In our analysis, we adopt the paradigm that lepton number violating neutrino-antineutrino 
couplings B are small compared to lepton number conserving terms H, H. We shall, in fact, 
adopt B as a small quantity for a perturbative solution [10j for the time evolution. To this end, 
we decouple the evolution of neutrinos and antineutrinos from the evolution of the neutrino- 
antineutrino system by introducing a (subsidiary) time evolution operator G(t) for the neutrino 
and antineutrino systems. The latter is then used to transform the time evolution equation for 
U(i, to) to what we shall henceforth refer to as the interaction picture. Establishing this line of 
action, we write 

V(t,to) = G(t,t )V l (t,to), (15) 

with 

^-G(t,t ) = -m (t)G(t,to), (16) 
at 
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where the subscript I indicates the associated quantity in the interaction picture. It is readily 
seen that the time evolution equation in the interaction picture is given by 

^Ui(Mo) = -ffl T (i)Ui(Mo) , (17) 



with 



Since the Hamiltonian Ho is lepton number conserving, the solution of Eq. (|16p is obviously of 
the form 



Hi(t) = G" 1 (t,to)'5H(i)G(t,to)- (18) 



Putting this solution into Eq. (I18D . we obtain 

M(t) ( , ° G-Ht,t )B(t)G(t,t )\ 

m) ~ {G-\t,t )BHt)G(t,t ) J' (20) 

We can now perturbatively solve the time evolution equation in the interaction picture by 
means of the Magnus expansion We write its solution as a matrix exponential 

U l (t,t ) = e n ^\ (21) 

where the Magnus operator Q(t,to) is the sum of the so-called Magnus approximants f2i(i,to) 
according to 

Q(t, t ) = Oi(t, t ) + 2 (t, t ) + . . . , (22) 

which are small in an appropriate sense as we shall see in due course. Various methods have 
been worked out to calculate the Magnus approximants [12]. They are found to obey 

Oi(Mo) = -i I drMi(r), (23) 
J to 

n 2 (Mo) = ~\ I d *i / d *2 K(ti),Hi(t2)l, (24) 
1 Jto J t 1 i 

up to second order in the perturbation. We shall see in our analysis that neutrino-antineutrino 

oscillation effects only enter the oscillation probabilities in second order perturbation theory. 

Our aim is to comment on how lepton number violating effects manifest themselves. We 

assume that the lepton number conserving effects, encoded in G, can be obtained either exactly 

or perturbatively by solving the underlying Schrodinger equation. The solution for the time 

evolution operator V(t, to) can then be written as 

U(Mo) = G(Mo)e c( *'' o) , (25) 

according to Eqs. (I15p and (I21D . Inserting the block form of the interaction picture Hamiltonian 
from Eq. (|20p into Eqs. (|23|) and (|24p. we find that the Magnus operator up to second order 
perturbation theory can be written as 

where we have introduced the shorthand notations 

,-t 



B(t,t ) = [ dr G -1 (t, to)B(t)G(t, to), (27) 
C(t,t ) = -\f t ^J^^ jfifa^Vz) - {b^B {t 2 )j^ , (28) 
D(t,t ) = ~j\t x j\tz Ib^B^)- ^S + (*i)S(t2)^ 1 ■ (29) 
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The dot indicates a derivative with respect to time. Note that the quantities C(t, to) and 
D(t,to) are quadratic in B(t,to), i.e., quadratic in the lepton number violating parameter B in 
the Hamiltonian. Then, up to second order terms in B, we obtain 



1 + 0( 2 ) (Mo) -*£(*, to) 
-w\t,t ) i + 6W(t,t ) 



+ ... 



where 



(2) (t,t o ) = C^-^Mo^Mo), 
(2) (Mo) = ^^-^(Mo^Mo)- 



(30) 
(31) 

(32) 
(33) 



These quantities show how lepton number violating neutrino-antineutrino mixing encoded in 
B(t) affects the neutrino and antineutrino sectors respectively. 

With these ingredients, we can now write down the oscillation probabilities of various kinds. 
Let us denote the oscillation probabilities in absence of lepton number violating terms by the 
notations 



P (u a ->■ v b ;t) 



Gbait) 

G ba (t) 



(34) 



Up to second order in lepton number violating terms, the probabilities would then read 

P(u a -> v b -t) = P Q (u a -> v b ;t) + 2Re ^G ba (t) x (c{t)0^ (t)) } , (35) 
P{y a -> u b ;t) = P {9 a -> z7 6 ;t) + 2Re{GL(t) x ^(^©^(t))^} . (36) 



Note that there is no summation over recurring indices implied here or elsewhere in the pa- 
per. Thus the lepton number preserving oscillation probabilities get modified. On the other 
hand, there can also be lepton number violating oscillations now, and the probabilities of such 
oscillations are given by 



P(v a -4 v b \ t) 
P(u a -4 u b ; t) 



G(t)B'(t) 
G(t)B(t) 



ba 



bu 



(37) 
(38) 



These oscillation probabilities for lepton number preserving as well as violating neutrino (an- 
tineutrino) oscillations Eqs. (135H38j) present the main result of our analysis. A few comments 
are in order. 

First, we want to emphasize that this method for obtaining the oscillation probabilities does 
not involve the effective mixing matrix at any stage of the analysis. It is of course possible, in 
principle at least, to diagonalize the entire Hamiltonian including the lepton number violating 
terms and derive oscillation probabilities from the resulting eigenvalues and eigenstates. Our aim 
was to demonstrate how some of the essential features of such a solution can be obtained without 
performing such a diagonalization. Note also that this approach reduces a 2N x 2N problem 
(./V being the number of neutrino species involved), the diagonalization of the Hamiltonian, to a 
N x N problem, since the perturbation expansion only involves block entries from the original 
Hamiltonian W(t). 

Truncating the exponential of the Magnus operator results in loss of unitarity. Since we 
truncate the series in a way that we only keep terms up to second order in the small perturbation 
B(t), we also expect that unitarity of the time evolution operator is only conserved up to second 
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order in this quantity. It is straightforward to verify this by invoking the unitarity condition on 
U(i, to). Note that such considerations also pertain to solutions for the neutrino and antineutrino 
evolution operators G(t,to) and G(t,to), if their solution is obtained in a similar way, i.e., by 
truncating the exponential of the associated Magnus operator in a perturbative manner. The 
fact that we seek a perturbative expansion of the solution to the time evolution operator readily 
implies that the structure of the oscillation probabilities is of the same kind. However, it should 
be realized that the probabilities given here are viable approximations for times that satisfy 
\B\t <C 1, where \B\ indicates the magnitude of any non-zero eigenvalue of B. Hence the loss 
of unitarity is small as it is triggered through the small quantity B(t), and will appear only 
at orders B 3 or higher which have been neglected in writing down the expression for various 
oscillation probabilities. 

In fact, unitarity can be tested explicitly from the expressions of the oscillation probabilities 
given above. Take, for example, Eq. (i35j) . It shows that 

E%^^*) = l + 2Re^{G^(i)x (G(t)0M(t))J 
b b a 

= l + 2Re(G t (t)G(t)0 (2) (t)) = l + 2Refo (2) (t)) , (39) 

V J aa V / aa 

using the unitarity of the matrix G. Further, looking at the expression of OW(t) in Eq. ([32]), we 
see that the diagonal elements of this matrix must have negative real part, since the diagonal 
elements of C are purely imaginary and those of BB* must be real and positive. Thus, start- 
ing from a particular flavor of neutrino u a , the total probability for lepton number conserving 
oscillations is 

Y^P{v a ^u h ;t) = l-(B{t)B ] {t)) , (40) 

L — ' V / aa 

b 

which is less than unity, as one should expect. It should also be noted, from Eq. (13T|) . that with 
the same initial state, the total probability for lepton number violating oscillations is given by 

= £ (SflSW (G(t)B\t)) = (B(t)G\t)G(t)B\t)) 
b b a a aa 

= (B(t)B\tj) , (41) 

V / aa 

using the unitarity of the matrix G. Therefore 

^2P(u a ^u b ;t) + ^2P(u a ^u b ;t) = l, (42) 
b b 

as expected, implying that in presence of lepton number violation, the lepton number conserving 
oscillation probabilities should decrease, making room for lepton number violating oscillations. 

The form of the oscillation probabilities given in Eqs. (I35H38P is illustrative in the sense that 
for the limiting case of vanishing neutrino-antineutrino coupling B(t) — > it readily reduces 
to the standard, lepton number conserving neutrino oscillation results. The appearance of 
the oscillation probabilities for v a — > v b and v a — > v b can also be interpreted in an intuitive 
way: reading Eqs. (1371 - 138]) from right to left they state that B (t) (B(t)) switches the initial 
neutrino (antineutrino) state to the associated antiparticle and G(t) (G(t)) then evolves the 
antiparticle (particle) state until its detection. So neutrino-antineutrino oscillations are clearly 
a signal for lepton number violation in the neutrino sector. Note that this statement is a general 
one: all oscillation probabilities come with B(t) and derivatives thereof; put another way, in 
an approach in which neutrino-antineutrino couplings are treated as a small perturbation to all 
other potential enhancements in the neutrino and antineutrino sector respectively, one cannot 
have modifications in the neutrino-neutrino and antineutrino-antineutrino probabilities without 
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introducing neutrino-antineutrino conversions at the same time. The oscillation probabilities 
for v a — > i>b and v a — > are generically different by virtue of the generic difference between 
G(t) and G(t). Even if the respective Hamiltonians H and H do not discriminate between 
particles and antiparticles (e.g. H(t) = H(t)) there still is a difference due to the fact that 
neutrino-antineutrino coupling can be complex and not self-adjoint, i.e. B(t) ^ B (t) in general. 

A similar situation pertains to u a —> and v a —> Uf, oscillations. They differ in the time 
evolution operators for neutrinos and antineutrinos, but also in the second order operators 
0^\t) and CK 2 ) (£). Those operators are, in general, not identical, which is again due to the 
fact that B(t) does not have to be self-adjoint. This statement translates to the fact that lepton 
number violating neutrino oscillations discriminate between particles and antiparticles even if 
there is no difference in the respective Hamiltonians for neutrinos and antineutrinos. This 
behavior is also expected from lepton number violation; note, that this distinguishing feature is 
generated dynamically rather than being imposed by hand. 

Having established an approach to neutrino oscillations in which lepton number violation is 
considered a small effect, let us now elaborate on the physical notion underlying this framework. 
For such purposes we begin by studying the one generation case of neutrino oscillations, since 
already in this case lepton number violation allows neutrino-antineutrino oscillations to develop. 
Making use of the formalism, we obtain 



B(t,t )= / drB(T)e iAi? M, (43) 

where we have defined 

AH(t) = ! dr AH(t) = f At [H(t) - H(r)] (44) 



as the difference between (CP non-conserving) potential terms in the neutrino and antineutrino 
sectors. Clearly, if the difference in such potential terms vanishes for some time t = t res , the 
integral Eq. (|43h has a stationary phase and we can evaluate it by means of the saddle point 
approximation. This yields 



B~e iAH ^ J- — . (45) 

V 2 7res 

Here we have introduced the adiabaticity parameter at resonance 7 res , which we shall properly 
define and explain shortly. To this end, let us start from the Hamiltonian of Eq. ([2]) and notice 
that for a simple two dimensional case, we can always easily find a unitary transformation (for 
a two dimensional problem phases are irrelevant and the unitary transformation amounts to a 
time-dependent rotation in flavor space), which entails an effective mixing angle Q(t) fixed via 

cose(f ) = MM sM , 1)= M, (46) 

"eft ^eff 

The effective oscillation frequency u e s of the system is here given by 



u; eff = ^4|B(t)| 2 + Atf 2 (t). (47) 

It is interesting to note that the effective mixing becomes maximal, neutrino conversions undergo 
a resonance, in the case of vanishing AH(t), i.e. AH(t rcs ) = at some resonance time t Tes . Put 
another way, mixing between neutrinos and antineutrinos becomes maximal when the difference 
between potential terms for the two species is minimal (or in fact vanishing at the resonance 
time). This is also corroborated by the fact that for identical neutrino and antineutrino potentials 
mixing is always maximal. Obviously, for the case of CP conserving neutrino and antineutrino 
potentials the difference between those vanishes at all times since they are identical to begin with. 
However, in the case of CP non-conserving matter potentials, the difference is generally non-zero 
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and can be varying with time also — the latter effect is crucial for some physics implications of 
the resonance structure as shall be seen in due course. It is also conceivable to allow for CPT 
violating potential terms, which then generate a non- vanishing difference AH(t). 

With the effective mixing at our disposal, we write the adiabaticity parameter at resonance 

as 



This definition of an adiabaticity parameter can also be understood physically as follows: the 
characteristic time of the neutrino- antineutrino system is r sys ~ l/cj e g, whereas the characteristic 
time of the interaction can be given as rjnt ~ (dO/di) -1 . Hence a small 7 res states that the 
system's time scale is much smaller than the interaction's time scale. Put another way, the 
transition is adiabatic. 

For 7 rcs > 1 we encounter non-adiabatic neutrino- antineutrino conversions; 7 rcs <C 1 gives 
the adiabatic case. So roughly speaking (neglecting the derivative of the potential term for 
the time being) the smaller B(t), the larger the non-adiabaticity of the neutrino-antineutrino 
system and hence if we start the evolution with only v a {pa) states present the transition u a — > u a 
{p a — > v a ) prevails; the lepton number violating oscillation channel v a — > v a (y a — > v a ) gets more 
and more suppressed as the non-adiabaticity increases. If, however, the change of the difference 
in matter potentials is small with time, this effect can partially compensate a small neutrino- 
antineutrino coupling, driving the evolution of the system towards adiabatic transitions opening 
the lepton number violating oscillation channel u a —> v a [p a — > v a ) again. Note also that the 
suppression of the lepton number violating oscillation channel depends on the time dependence 
of the potential terms in the Hamiltonian. 

Substituting the result for the adiabaticity parameter 7 res of Eq. (138]) in the expression for B 
of Eq. (|45p and at the same time keeping in mind that the oscillation probability P{v a v a \t) 
of Eq. (|37p is directly proportional to B , it is seen that interpreting lepton number violating 
(neutrino-antineutrino) couplings as a small perturbation to lepton number conserving neutrino 
oscillations in the one generation case is equivalent to assuming that neutrino-antineutrino os- 
cillations occur non-adiabatically; non-adiabaticity hence closes the lepton number violating 
oscillation channel v a — > i> a , but at the same time improves the perturbative expansion as 
outlined above. 

These results for the one generation framework make the case for referring to the perturbation 
theory as developed in this letter as a non-adiabatic perturbation expansion. 

From the definition of the adiabaticity parameter 7 res in Eq. (148 p another interesting feature 
emerges. Suppose we assume small lepton number violating couplings B(t) between neutrinos 
and antineutrinos in the oscillation Hamiltonian. This means that oscillations between particles 
and antiparticles are suppressed by the large non-adiabaticity of the transitions giving rise mostly 
to the lepton number conserving oscillation channel u a — > v a . If, however, the time variation of 
the difference in matter potentials at the resonance is sufficiently mild, it is seen that the presence 
of such matter along the neutrino (antineutrino) propagation path can drive the system towards 
adiabaticity thus opening the oscillation channel v a — > v a . Put another way, the presence of 
CP (or even CPT) non-conserving matter of a varying density (readily re-interpreted as a time 
dependence of the potential terms) can enhance lepton number violating neutrino oscillations 
as compared to the case in vacuo. Note that this statement holds regardless of the nature of 
the perturbation theory developed in this letter, since it is merely a result obtained from the 
adiabaticity parameter 7 res and the latter does only need the Hamiltonian of the system as a 
prerequisite. 

Let us next approach the two generation case in vacuo. In this case, apart from a term 
proportional to the unit matrix, the flavor oscillation of neutrinos as well as antineutrinos is 




1 d0(t) 



1 



dAH(t) 



(48) 
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governed by the Hamiltonian 

Am 2 



H=H 



4p 



o x sin 26 — a z cos 26) = uj ( a x sin 26 — a z cos 26 



(49) 



where the it's are the usual Pauli matrices, and 6 is the mixing angle in absence of lepton number 
violating terms. Accordingly, the evolution operator in the neutrino as well as the antineutrino 
sector will be given by 

G(t,t = 0) = G(t,t = 0) = cos ujt — i (a x sin 26 — a z cos 29^ sin ujt , (50) 

in absence of lepton number violating terms. In order to include the effects of the lepton number 
violating terms, let us first write the matrix B(t) as 



B(t) = bo(t) + bi(t)ai , 



where 



bo(t) = -trB(t) 



&i(t) = -tx[aiB(t) 



(51) 



(52) 



Note that if CPT is conserved, Eq. ([7]) dictates that we must have b y (t) = 0. It is straightforward 
to see that 



where 



A>(*) 



B{t) = /3 (t) + Pi(t)<Ti 



bo(t), 

sin 26 b x (t) sin 26 -b z (t) cos 26 + cos 2ut cos 26 
—b y {t) sin 2ut cos 26 , 

b y (t) cos 2ut + sin 2cot b x {t) cos 26 + b z (t) sin 26 



(53) 



b x {t) cos 26 + b z (t) sin26 



cos 26 



b x (t) s'm26 -b z (t) cos 26 



+ cos 2u;t sin 20 



b x (t) cos 26* + 6 z (i)sin26 l 



-b y (t) sin 2cjt sin 20 . 



(54) 



We will assume here that the elements of B are time- independent. These expressions can then 
be easily integrated to give B, and the resulting expressions plugged into the relevant formulas 
to obtain various oscillation probabilities. 

For the most general B, the expressions are cumbrous. Here, we present an illustrative 
example for a very simple form of B which can be analytically tackled without much trouble. 
Consider B to be the multiple of a unit matrix, which means that the lepton number violating 
effects are somehow flavor blind. In Eq. (|5ip . this would mean that all i>, = 0, only bo is non-zero. 
Now, from Eq. (|54p . we find that B is b^t times the unit matrix. Looking back at Eqs. (|37p and 
(I38|) . we find that the probabilities of neutrino-antineutrino oscillations are given by 

(55) 



P(y a -> u b ;t) = \b \ 2 t 2 x P (v a v b ;t) , 
whereas the lepton number conserving oscillation probabilities are given by 
P(u a -> u b ;t) =(l- \b \ 2 t 2 ) x P {u a -> u b -t) . 



(56) 

Let us discuss the nature of this solution. If there is no lepton number violation, a plot 
of oscillation probability vs. t would exhibit periodic highs and lows of the same magnitude. 
If there is lepton number violation, these periodic nature would be modulated by the function 
(1 — |foo| 2 i 2 ). This means that, as the path length increases, the successive maxima and minima 
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of probability would be less and less pronounced. There is, of course, no fear of the maxima 
and minima vanishing and the curve becoming completely flat, because, as said earlier, these 
solutions should be valid in the region where \bo\t <C 1. 

In fact, with this particular choice of the matrix B, one can do better and solve the problem 
exactly. To see this, let us recall the form of the interaction Hamiltonian given in Eq. (|18p . 
Since G = G for the Hamiltonian Mo, we obtain the 2x2 block structure of Hi as follows: 

o)=N( e -i (57) 

where bo = |6o|e JQ . Moreover, in the exponent in the Magnus expansion, only the term Oi is 
non-zero. From Eq. (j23[) . we find that it is given by 

n(t) = «i(t) = -ffi.it. (58) 

Because of the special form for the matrix Hj, the exponentiation can be done exactly, and one 
obtains 

e°W = (j j)cos|6 |t-*( e iQ e o a )sin|6 |t. (59) 

The evolution operator of the system is defined through Eqs. (|15p and (|2ip . It shows that 
P{y a -)• v h \ t) = cos 2 \bo\t x P (u a -> z/ 6 ; t) , 

P{v a ^v b -t) = sm 2 \b \txP (v a ^i; b ;t). (60) 

The expressions given in Eqs. (|55p and (j56[) are nothing but the approximations up to leading 
order terms in the time dependence. However, in the form given here in Eq. (I60p . it is valid for 
all t. Note that, although this last illustration uses a time-independent Hamiltonian, the general 
formulas that we give in Eqs. (1351 138p are valid for time-dependent cases as well, e.g., in matter 
induced oscillations involving neutrinos and antineutrinos. 

In conclusion, we have developed a practical and efficient way of dealing with lepton number 
violating neutrino (antineutrino) oscillations. In a non-adiabatic perturbation theory, in which 
lepton number violation is treated as a small effect compared to common flavor oscillations, 
as is compatible with experimental findings, we have given explicit expressions for the vari- 
ous oscillation probabilities up to second order in the perturbation and have embarked upon 
their interpretation using illustrative examples for one and two neutrino generations. For the 
case of one neutrino generation, we have found a novel resonance in lepton number violating 
oscillations, which in CP non-conserving matter enhances a potentially small lepton number 
violation in vacuo. For the exactly solvable two generation case of time-independent and flavor 
blind lepton number violation in vacuo, we understand that in the presence of lepton number 
non-conservation the common flavor oscillation probabilities are periodically modulated with a 
frequency given by the lepton number violating coefficient. 
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